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SUMMARY 


Correlation studies show that some degree of coherence exists through- 
out the whole of the observed plasma turbulence spectrum, although time 
and length scales of coherence decrease markedly with increasing frequency. 
Propagation characteristics and energy dissipation mechanisms of these 
turbulent eddies are being investigated. We find that growth of a strongly 
coherent drift wave instability is associated with an increase in amplitude 
of the high frequency end of the turbulence spectrum. On the other hand 
growth of another type of coherent oscillation (an ion acoustic wave) reduces 
the amplitude of the turbulent spectrum, but only in the frequency regions 
adjacent to that of the strong coherent oscillation. Further experimental 

j 

data is required before these results can be related to theoretical work. 



I. Introduction 


Turbulent fluctuations of electrostatic fields and charged particle 

18 

density in a plasma medium are connected with plasma transport properties. 
However the study of plasma turbulence is even more complex than the study 
of fluid turbulence, because charged particles give rise to collective 
oscillations in a plasma medium. Collective oscillations are rarely observed 
in a fluid. 

By inserting a Langmuir probe in the plasma and displaying the resulting 
signal on a spectrum analyzer we observe the oscillations in plasma poten- 
tial if the probe is at floating potential, or the oscillations in charged 
particle density if the probe is biased into the ion saturation region. 

All data presented in this, report were obtained by use of probes at floating 
potential. 

Our spectrum analyzer can be used to display the amplitude of the 
signal versus frequency for frequencies up to 1 MHz. The spectrum is 
continuous over this range, see Fig. 1-1. By varying some of the external 
parameters associated with the operation of the arc such as background 
pressure, gas flow through anode and cathode, magnetic field strength and 
magnetic field geometry, we can cause the onset of large amplitude coherent 
oscillations, see Fig. 1-2, and also affect the level of the turbulent 
oscillations. Although the rest of the oscillations are lower in amplitude, 
they occupy a larger portion of the spectrum than the coherent oscillation. 

We are, therefore, interested in the effect of the onset and growth of these 
instabilities on the complete spectrum of oscillations. 
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Calculation of an auto correlation function for specific oscillation^ 

leads to a power spectrum. By monitoring changes in the power spectrum 

associated with changing external plasma parameters we can establish the 

form of coupling between the large coherent instabilities and the turbulent 

portion of the spectrum. It has already been shown that anomalous plasma 

diffusion across the magnetic field is related to the presence of large 

1 2 ? 6 12 

coherent instabilities 5 ’ 5 ’ , and the study of these coupling effects 

should give useful information about the energy exchange process involved 
in the plasma. 

Using this correlation technique we were able to detect some degree 
of coherence for oscillations throughout the frequency spectrum up to 
100 KHz, and also obtain order of magnitude estimates of the temporal 
extent of these oscillations. Calculation of cross correlation functions 
for different oscillations in the plasma was used to obtain propagation 
characteristics of these oscillations and the order of magnitude of the 
spatial extent of these fluctuations. 

II. Power Spectra 

The recent acquisition of new instrumentation allowed us to measure 
4 l4 15 

power spectra ’ ’ , that is, to plot power contained at each oscillation 

frequency of the plasma, in a more reliable way than can be done from a spec- 
trum analyzer. The PAR Model 102 Fourier Analyzer gives the power spectrum 
directly from an auto correlation function of the fluctuations detected by 
a Langmuir probe. This instrument is designed to make use of the Weiner 
autocorrelation theorem, which states that the Fourier cosine transform 
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of the auto correlation function C^( r ) gives directly the power spectrum 
of oscillations, 0^(<u) : 

oo 

011 (u>) =2 G^( r ) cos COT dr 2-1 

o 

When external plasma parameters such as magnetic field, neutral gas 
pressure, and gas flow rates through both anode and cathode are varied, 
the observed power spectra exhibit definite changes. By establishing the 
connection between the characteristics of the turbulence spectra, the 
development of strong coherent oscillations, and the internal parameters 
of the arc, we can therefore investigate energy exchange mechanisms experi- 
mentally. 

Theoretical work along these lines has been carried out, for example, 
by Tchen^ 5 He has considered the spectral distribution of energy for 
charged-particle density fluctuations in a plasma in a strong magnetic 
field, for both the collisionless and the collision-dominated case. 

Describing the plasma by the Wavier -Stokes equation of motion including 
an electrostatic field, the equation of continuity for a compressible 
fluid, and the Poisson equation for charge densities, he solved the 
resulting non-linear equation using an approximation called cascade decom- 
position. The Fourier components of all eddies were grouped into two 
separate parts: large-scale eddies, which contribute to the development 

of the energy spectra, and smaller eddies which set up turbulent properties. 
The turbulent motion of the big eddies was considered homogeneous, stationary 
isotropic and compressible. The small eddies were assumed incompressible, 
non- homogeneous, anisotropic, non- stationary and did not contain as much 
energy as the large scale eddies. 
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Two different energy exchange mechanisms had to be assumed to account 
for the modal transfer from big eddies to smaller ones. Dissipation, or a 
drain of energy was assumed to occur for the collision dominated plasma by 
a molecular viscosity, and for the collisionless plasma, by collective 
electrostatic fluctuations. The spectra derived under these assumptions 
depend on the wave number k of the oscillations raised to an appropriate 
power. 

For the colli sional plasma the density fluctuation spectrum can be 
represented as 

G OCk -3 / 2 k < kg 

or 

GoCk -9 / 2 k > kp 

where k^ = (4 a/a X 2 )^ 2 
a = (e ay^) 1 / 2 
e = a dissipation quantity 
X = coefficient of Bohm diffusion 
<x> c = cyclotron frequency 

For the collisionless plasma the density fluctuation spectrum 

G<OCk“ 5 

These expressions show how the type of energy dissipative mechanisms affect 
the shape of the spectrum of density fluctuations. 

In order to experimentally test the results of Tchen's theory it is 
necessary to isolate various wave- number components and find the power 
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present in eacb component. However techniques for the isolation of parti- 
cular wave number components are not available at this time, and in the 
present experiment particular frequency intervals were isolated. Therefore 
a method of converting from a frequency to a wave number spectrum is needed 
in order to test theoretical expressions which have been derived. One possible 
technique for converting from frequencies to wave number (to — ►k) is proposed 
in a later section. 

In our experiment three different operating regions of the plasma were 
selected and the effect of neutral gas pressure on the power spectrum was 
studied: 

(i) a quiescent regime of the plasma with neutral gas pressure 

-4 

p = 1.9 x 10 torr, main solenoid magnetic field B = 770 gauss, and the 

magnetic field in the neighborhood of the cathode, B = 3^5 gauss. 

c 

(ii) a regime where, by reducing the flow-rate of argon through 
the cathode 20%, onset and growth of one coherent instability, an ion 
acoustic wave, at about 4 KHz in the frequency spectrum, was initiated. 

(iii) a regime where the onset and growth of another type of 
coherent instability, a drift wave, at about 7 KHz in the frequency spec- 
trum was initiated. This was accomplished by increasing B to 2100 gauss 
and adjusting B q to the same value. 

The effect of pressure change in region (i) is shown in Fig. 2-1, The 

-4 

neutral gas pressure was increased from an initial value of 1.9 x 10 torr, 
first 20$, and then 6o$, by closing the baffles between the system and the 
diffusion pumps. The following changes in the power spectrum are observed: 
at the low frequency end, magnitudes are affected by up to 100$ but there 




Figure 2-1 Power Versus Frequency ^ 

Power in arbitrary units - initial pressure P = 1.9 x 10 torr. 

b h = 770 s auss ’ B c ,= 385 gauss 
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is no systematic change with increasing pressure. At the high frequency •* 
end the slope of the power spectrum is unchanged. 

The effect of pressure change in regime (ii) is shown in Fig. 2-2. 

It should be noted that initially the power contained in the coherent 
oscillation at around k KHz is at least 60 times the power contained in 
adjacent frequency components of the turbulent spectrum. By increasing 
the pressure in the same way as before, the amplitude of the ion acoustic 
wave can be reduced to the point where the wave is no longer distinguish- 
able from those oscillations around it. The two power spectra (with and 
without the ion acoustic wave) are shown in Fig. 2-2 and it is clear that 
the power in the oscillations at .frequencies on either side of the ion 
acoustic wave increases as the power in the coherent oscillation decreases. 
However the high frequency end of the spectrum is unaffected. This suggests 
that the ion acoustic wave gains some of its energy at the expense of those 
low frequency oscillations which surround it on the frequency spectrum. 

The effect of pressure change in regime (iii) is shown in Fig. 2-3. 

In this case the strong oscillation at around 7 KHz is a drift wave, and 
the power contained in this oscillation is about 60 times the power con- 
tained in adjacent frequencies of the turbulent spectrum. It can be seen 
in Figure 2-3 that reducing the amplitude of this coherent instability by 
increasing the pressure by 6.0$ has little effect upon the low frequency end 
of the power spectrum. However there is a noticeable change in the slope 
of the power spectrum above 50 KHz. 

The slope of power spectra 2-1 through 2-3 are all different, which 
suggests that different energy exchange mechanisms are present when a 




Figure 2-2 Power Versus Frequency. (ion acoustic, wave) Power in arbitrary 

units - initial pressure P = 1.9 x 10 , = 770 gauss B c = 385 gauss 



Figure 2-3 


Power Versus Frequency. (dri^t^wlve) Power in arbitrary units - 
initial pressure P = 1.9 x .10"^ torr, B H = 2100 gauss j B c = 2100 gauss 
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large coherent oscillation is present in the spectrum of oscillations, and 
that the type of instability determines which process is dominant. 

Our results show therefore that (a) both the amplitude and the shape 
of the power spectrum depend on the operating conditions of the arc. (b) 

The shape of the power spectrum depends on whether or not an instability 
is present, (c) The three instabilities investigated must involve different 
energy exchange mechanisms. However no clear criterion is available for 
distinguishing between collisional and collisionless regimes in the plasma. 
Therefore unambiguous interpretation of our results is difficult at this 
stage of the investigation. 

III. Degree of Coherence in Turbulent Portion of the Spectrum 

The term" coherent oscillation" in the following discussion means an 
oscillation which can be detected over a finite spatial region of the 
plasma, and is observable for a finite time. This is not a precise defini- 
tion, but we use it to distinguish between those charged-particle density 
and electric-field fluctuations which occur completely at random in the 
plasma, and those in which a definite phase relation exists between syste- 
matic perturbations in these quantities, measurable over some definite 
time interval and some observable region of the plasma. 

On a spectrum analyzer the observed Langmuir probe signals appear to 

represent random noise from the plasma, but our results show that some degree 

■£ 

of coherence does exist throughout the frequency spectrum up to 100 KHz . 

This is the upper limit for using our correlation technique with the 

present instrumentation. 

>■ 
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An auto correlation function compares a signal with itself over a 

length of time T . Mathematically this is expressed as 

T 

C i;l ( 7 ) = 1/2T / f 1 (i ) f x ( t + r ) dt 3-1 

-T 

For example the auto correlation function of wide band noise (a purely 
random signal with no coherence) can be seen in Figure 3~1» The magnitude 
of the auto correlation function at “f = 0 on the curve is a measure 
of the magnitude of the mean square value of the amplitude of the oscilla- 
tions in the wide band noise 

T 

c n(°) = 1/2T J ) dr 3-2 

-T 


For times greater than - 0, the auto correlation function falls to zero 
This is due to the fact that signal f^(*t ) is purely random, with no coher- 
ence. Therefore after f^("t) is displaced an amount T , f^(*fc}_ and 
f^( 'fc' + T ) are no longer coherent with one another. 

On the other hand the auto correlation function of a pure sinusoidal 
oscillation is itself a repeating sinusoid (Figure 3 -2 ). Mathematically 
this can be written in the form 



T 



mLirn 
— ► o>0 

1/2T f 

-T 

sin cd *fc. 

sin cd (-fc +T ) dt 

3-2 

Lim T- 

— * 

T 

T 

[ 1/2T l 

s i n cd ot- 

sin cDt: cos cdI* dt 

*1 

1/2T J 
-T 

>■ 

sin co T 

COS CD it 

sin cd *Y d ”t I 


= l/2 cos CD T* 
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An auto correlation function taken over a wide band of frequencies from 
the plasma has its greatest contribution from any coherent oscillation 
rather than the incoherent portion of the spectrum. For example, if the 
spectrum of oscillations consisted of wide band noise plus one coherent 
oscillation at a discrete frequency f , for times later than T = 0, C^( T' ) 
would become a repeating sinusoid at the same frequency f, and for large 
enough Y , the repeating sinusoid would decay as the discrete oscilla- 
tion began to lose coherence with itself. 

When the operating conditions of the arc are set so that no large 
instability at any particular frequency is present in the observed spectrum 
analyzer display, the auto correlation function of the resulting signal 
does not have the shape associated with wide band noise. Such an auto 
correlation function is shown in Figure ’3-3* This case is intermediate 
between the auto correlation function computed for wide band noise and 
the one for a pure repeating sinusoid. A possible interpretation of the 
form of ) is that many short term coherent oscillations at different 

frequencies are superimposed on one another at the spatial location in the 
plasma where the probe is located. We therefore instituted a new technique 
of restricting the bandwidth of the incoming signal from the probe. The 
signal from the Langmuir probe was passed through a narrow band-pass filter, 
adjusted to center on the frequency f to be studied. Since the signal 
level was low, a low noise wide-band amplifier was then used to increase the 
amplitude of the transmitted signal. This signal was then fed into a variable 
Q, selective amplifier, also centered on the frequency f, in order to obtain 
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a high enough amplitude for it. to be processed by the correlator. The 
final signal was then fed into the correlator which printed out the auto 
correlation function vs T (quasi-time) on an X-Y recorder. 

We found that even when the center frequency of the filter was set on 
a portion of the turbulent spectrum that was considered to be totally ran- 
dom in nature, the resulting auto correlation function showed a 

definite sinusoidal behavior, with amplitude decaying in time. This can 
be seen in Figure 3-*+« 

We must allow for the possibility that the form of this auto correla- 
tion function could be produced by instrumentation, e.g., the response of 
"ringing" filters, rather than by the plasma. However there is experimental 
evidence which shows that the plasma does determine the observed form of 
the correlogram: this involves the time, scales of coherence observed, and 

the phase shifts which are measured. 

g 

The time scale of coherence, Y , is the time it takes for an 

c 

auto correlation function to decay to l/e of its value at Y = 0. If the 
observed auto correlation functions were due entirely to instrumentation, 
their relative magnitudes might be affected by plasma conditions in the 
region of the probe from which the signal is derived, but the time scale 
for decay of the correlogram should be independent of the operating condi- 
tions of the arc. Our results clearly show that the time scales are affected 
by internal plasma parameters. Figure 3-5 shows a plot over one order of 
magnitude of time scales of coherence versus frequency, measured from the 
correlograms obtained over the plasma turbulence spectrum, and plotted for 




Figure 3“ 3 Auto correlation Function of Signal detected by a Langmuir Probe 

at Floating Potential placed in the Plasma, (No large instability 
present). 



Figure 3-^ Auto correlation Function of Signal detected by a Langmuir Probe 

at Floating Potential placed in Plasma (Large instability present). 
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rnany different operating regions of the arc. It can be seen by comparing 
these points that X varies as much as 25 $ under different arc condi- 
tions. Even if part of the auto correlation function computed is due to 
the instrumentation involved, a measure of the relative change in the time 
scale of coherence at each frequency due to the plasma itself can be calcu- 
lated from C^( T) . We are developing a computer program to separate out 
the true coherence decay time scale from the correlograms in order to 
obtain reliable absolute values for X toe to the plasma alone. 

Although it is not clear whether some of the calculated decay time is 
associated with the process of filtering the signal, the consistent measure- 
ment of changes in phase at different spatial positions shows that the 
plasma is determining the form of the correlogram. This point is discussed 
more fully in a later section. 

The length of time an oscillation remains coherent with itself is 
directly connected to energy dissipative mechanisms taking place in the 
plasma and the time scale for an oscillation to be modified or completely 
destroyed, shows how long energy is remaining in an oscillation at a 
specific frequency. We therefore investigated how the coherence time of 
oscillations in the turbulent .portion of the frequency spectrum depended 
on the presence of a large instability. Figure 3-6 is a plot of coherence 
time of oscillations versus frequency where two separate cases are shown, 
one with an ion acoustic wave present in the spectrum of oscillations, 
the other with this instability no longer' observable on the spectrum 
analyzer. Clearly the coherence times involved do not depend on the 
presence of the ion acoustic wave; there is no coupling, for example, between 
the ion acoustic wave at 4 KHz and the oscillations between 15 KHz and 100 KHz. 
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Figure 3-6 Time Scale of Coherence versus Frequency for two Plasma arc 
operating conditions (with and without an ion acoustic wave 
present-; 
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In summary, the quantity T" c , which is a measure of how long an 
oscillation survives, shows that (a) there is a degree of coherence in 
some oscillations in the frequency spectrum up to at least 100 KHz. (b) 
values of *T over the turbulent spectrum are independent of the presence 
of a strongly coherent oscillation. 

A further point to be considered is that when our probe is located 
at one position in the are, the observed decaying auto correlation func- 
tion may be due to two separate plasma effects. The first is that the 
oscillation is being modified as time goes by and is losing coherence with 
itself. The second factor affecting the apparent decay rate is convection 
of the oscillation past the single probe. If an oscillation in the plasma 
is moving past the stationary probe with some fixed velocity, with different 
amplitude in different regions of the plasma, any apparent changes detected 
by the probe will 'not be truly representative of the local region of the 
plasma in which the oscillation occurs. We need therefore to obtain a 
"moving-frame auto correlation function", which is a calculated auto-corre- 
lation function that simulates an auto correlation function taken from a 
probe moving at the convection velocity of the oscillation. In order to 

find the moving frame auto correlation function we have to measure cross 

7 

correlation functions of a number of spatially separated probes. This 
method of measuring X” c for an oscillation that is moving past a stationary 
probe is discussed in a later section. At the present time our measurements 
do not allow for this effect. 

IV. Propagation Characteristics 

The important propagation characteristics of plasma oscillations, such 
as wave number, wavelength, and speed of propagation (phase velocity) , may 
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be determined by monitoring the cross-correlation functions of signals 
detected by two spatially separated Langmuir probes. The cross-corre- 
lation function is obtained by passing the detected signals into a PAR 
Model 100 Signal Correlator which performs electronically a close approxi- 
mation to the mathematical definition: 

T 

C 12 ( r ) = i J f x (t) f 2 ( i- + r ) dt 4-1 

-T 

Therefore the cross correlation function compares two separate signals over 
a length of time . This is not the same as for the auto correlation 
function which compares a signal with itself. 

If both f^( *6) and f 2 (lfc) are purely random wide band noise, the 
resulting cross correlation function will be zero for all times since 
there will be no correlation between f^(-fc) and ,fg (■£:)., even at r = o, 
because the integral 

T 

Ci 2 (°) = lim T CO l/2T f f^i) f 2 (i) dt 4-2 

-T 

does not represent a mean square value of either function. 

For the case of two pure . repeating sinusoids of the same frequency 

the cross correlation function will be a repeating sinusoid. The value 

of C^ 2 ( X ) at T = 0 is determined by the phase relation between the two 

sine waves. If the two waves are represented by sin (cot + 0^ ) and 

sin (cot + 0 2 ) where 0^ and 0 2 are the phase angles involved, it can be 

shown that the cross correlation function 

T 

C AT) = Lim l/2T f sin (cot + 0 ) sin (cot + coT* + 0 p ) dt 4-3 

T —>£& "■ *m 
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can be written as 

C^g(T) = (l/2) cos coT cos 0 + (l/2) sin cd T" 1 sin 0 4-4 

where j t is the phase difference between the two sinusoids 

0 = 0 2 - 4 -5 

From 4-4 we can see that the value of C^g(T ) when T = 0 depends on 
the phase difference 0, 

C 12 (°) = S£|i 4-6 

Therefore if the two sine waves are in phase, (0 = 0 2 , 0=0), C^ 2 (o) will 
have its maximum value of l/2 at T = 0. If 0 = 0 2 - 0^ = 90°> C 12 (0) = 

and if 0 - l80°, C^g(O) = - l/2. Some graphical examples of this can be 

seen in Figure 4-1. 

If the two sine waves have a different frequency the orthogonality 

of the sine and cosine functions causes the cross correlation function to 

be zero for all time , since 

T 

C.p(T) = lim l/2T f sin cd. sin cnA'fc + T) dt = 0, cu. 0 CDp 4-7 

The two spatially separated Langmuir probes in the arc pick up signals 

and f 2 (’fe ), both of which are . fed into the signal correlator. If 
the same oscillation is being detected at both probes, the observed form 
of the cross correlation function will show a sinusoidal variation with 
time. if several different oscillations occur at both probes the 
resulting form of the cross correlation function will represent the super- 
position of the sinusoids that would result from each separate oscillation. 




Figure 4-1 Cross correlation function of f^-6) and fgC’t' ) for two separate 
phase differences. 
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Therefore a study of the cross correlation function of signals from 
two floating Langmuir probes gives information about properties of the 
electrostatic fluctuations in the plasma: 

1. If the cross correlation function displays a sinusoidal variation 
with T , the signal from probe 1 contains the same oscillation as the 
signal from probe 2. Therefore the size or spatial extent of this oscilla- 
tion must be at least the spatial separation of the probes. 

2. If a sinusoidal variation with T' is observed, the value of the 
cross correlation function at Y - 0 can be used to calculate .the phase 
difference ft in the oscillation at points 1 and 2. Further, if the spatial 
separation of the two probes is systematically changed, by monitoring the 
change in 0 one can calculate some of the propagation characteristics (wave 
number, wavelength, speed of propagation) of the oscillation. 

An illustration of the change in 0 with changing probe separation is 
given in Figure 4-2 for the case of a coherent instability of large ampli- 

g 

tude . Four different correlograms are shown, taken using two radially 
movable probes placed in the arc in a region of the plasma where a strong 
ion acoustic wave was observed on the spectrum analyzer. When these probes 
were aligned with one another ’ (curve No. l) the resulting correlogram had 
its maximum value at T” = 0 where the phase of the ion acoustic wave at 
both probes is approximately the same. Curves 2-4 shows that as one probe 
was moved radially away from the other probe (3 steps of l/2 cm displace- 
ment each) the value of C^( T) ’T = 0 is no longer a maximum. The 
maximum in the correlogram has shifted to the right as the probes are 
separated. This shows that the ion acoustic wave is propagating in the 
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radial direction. Determination of its phase velocity from measurement 
of radial phase shifts show a radial phase velocity of approximately 

. 4 / 

7.4 x 10 cm/ sec. 

The phase shift, A 0, with changing probe separation was measured 
by observing the first zero crossing of successive correlograms. From 
eqn. 4-4 at the point T* where the first zero crossing occurs 

= 0 = l/2 cos 00 cos 0 + 1/2 sin o> sin 0 4-8 

As the probe separation changes, the point at which the first zero 
crossing occurs at 2 / "j*^. The corresponding quantity which must 

also change in equation 4-8 is the phase difference 0. A simple graphical 
technique for determining the shift of the first zero crossing from 
V to 'T' g was used to calculate A 0. 

In order to see whether or not the same oscillation occurs at two 
spatially separated points in the plasma and to be able to calculate the 
phase difference, 0, of this oscillation at these two points, only one 
correlogram is necessary (see eqn. 4-4) . But if one wants to monitor a 
phase shift, A 0, with changing probe separation as indicated above, 
several correlograms are needed. The phase shift A 0 can thus be used 
to calculate propagation characteristics of the oscillation. Assuming the 
relationship 

0 = kl 4-9 

where k is the wave number of the oscillation, and 1 is the separation of 
the probes corresponding to a phase difference 0, the wave number k can be 
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directly calculated from the measured value of /. Note that three separate 
measurements are required for full information, since k^, k^ and k^ must 
he determined independently. 

If we know the wave number for separate frequencies in the spectrum 
(cd = co(k) ) we can plot the dispersion relation for the plasma point by 
point. From this curve both the phase velocity and the group velocity can 
be calculated 


Vp = co/k ; Vg = doo/ dk 
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In order to investigate the propagation characteristics of oscillation 

at selected frequency intervals throughout the turbulent portion of the 

frequency spectrum, the signals from each of the two probes was fed through 

separate identical filter .and amplifier circuits (see Fig. 4-3’). Just as 

for the case of the ion acoustic wave, when the radial separation between 

two probes is altered, the correlograms showed a phase shift, A 

In figure 4-4, which involves an oscillation at 15 KHz, as one probe 

is moved in relation to the other, the time 1 at which the maximum value 

of C^C'T’) occurs changes, and the first zero crossing of C^C^) shifts to 

the right. This indicates a phase shift A from which k^, and v^ for 

this oscillation can be calculated. It was found that A0 = 48° for l/2 cm 

displacement k = 1.18 cm \ \ = 3*74 cm and v = 5.6 x 10^ cm/sec. A 

X* P 

similar set of curves can be seen in figure 4-5 for an oscillation at 6 0 KHz. 
In this case a l/4 cm displacement led to a 50° phase shift k^ = 3*4 cm - '*', 

X = 1.8 cm and v = 1.08 x 10^ cm/ sec . Therefore for the ion acoustic 
wave at 4 KHz, and the oscillations at 15 KHz and 60 KHz our preliminary 
measurements show that in the radial direction the phase velocity decreases 
with increasing frequency of the oscillation. 
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Preliminary measurements have shown that consistent changes in the 
phase shift A 0 with radial separation only occur over small radial 
lengths. Therefore the calculated values for k^ depend on the radial 
positions at which the two probes have been placed. We interpret this as 
showing that only small radial sections or layers of the plasma can be 
considered homogeneous enough to yield interpretable results. Because of 
this, the dispersion relation cn = cn(k) is definitely not a simple relation- 
ship. 

V. Frozen and Non-Frozen Flow 

Care must be taken when interpreting correlation measurements in the 
presence of appreciable turbulence intensity. Shear forces, for example, 
may be continually creating and modifying the turbulence. In this case 
Taylor's hypothesis that turbulence may be regarded as a frozen pattern of 
eddies moving past an observer, is no longer valid. It has been shown for 
the case of fluid turbulence that Taylor's frozen flow hypothesis is valid 
only if the level of turbulence is low, viscous forces are negligible, and 
the mean shear is small"^. If one cannot assume a frozen flow, interpreta- 
tion of correlation measurements become more difficult. 

For the case of frozen flow at some convection velocity, a signal 
registered by a stationary probe at one point will be identical with that 
registered further "downstream" at- a second position x^ but will occur at 
a time ^ later. The cross correlation function G - | p (T ) will be a maxi- 
mum at *y = 7% and the velocity of convection, U , of the pattern may 

be defined as x^/ 7”" 




The cross correlation function can be written as 
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1 > ■ 


C 12 (x 


u 


t^) 
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If a number of cross-correlation measurements are performed at 

separations of x. ( i = 1, 2, 3> • • • • ) s each curve will pass through its 

maximum value when = T*. = x./U • The auto-correlation in a frame of 

1 v c 

reference which is moving with the convection velocity will be a straight 
line, always equal to a constant, namely the maximum value of the cross 
correlation function. The Weiner Auto-Correlation Theorem will lead to a 
power spectrum in which all the power can be found in an infinitely narrow 
band of frequencies centered at zero frequency. If the observer were to 
remain stationary, .the spectrum observed would be due to eddies moving 
past at the convection velocity and not due to any temporal turbulence 

7 

fluctuations . 

Any fixed- frame frequency f is related to the spatial extent, X , 
of the eddy producing the observed signal, by a relation of the form 


f * u 0 A 
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or 


cd = 2rtf = kU 

c 
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where k is the wave number of the eddy. Therefore the convection velocity, 
spatial extent, and wave number of eddies in a frozen flow may be determined 
from a knowledge of the spectrum of oscillations and a series of cross- 
correlation measurements. 

For the case of a non-frozen flow a series of cross-correlation curves 
will still rise to a maximum at some value of 'V , but this maximum value 
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will now be a function of x, the separation. This is because the flow is 
now being modified as it moves downstream. The convection velocity' can no 
longer be defined by the time delay at which the maximum of a particular ' 
measuring point separation. If an observer were to move at the convection 
velocity while measuring the auto correlation function, C_^( TO would, be 
the envelope of the cross correlation curves. The convection velocity 
would be defined by the time delay at which this envelope of all the cross 
correlation curves intersects the curve for a particular measuring point 

j 

separation (see Fig. 5-1). 

If an observer were stationary the auto correlation function he would 
measure would have two effects associated with it, l) The eddy is being 
modified 2) the eddy is being swept past. When discussing a time scale 
of coherence, only the modification of turbulence with time should be taken 
into account. Therefore for both the frozen and non-frozen pattern of 
turbulence the time scale of coherence as discussed in an earlier section, 
should be derived from a moving frame auto correlation. Note that in the 
frozen flow situation this auto correlation function will lead to a I 

c 

which is infinite, since the definition of a frozen flow is that the pattern 
of turbulence is not changing - with time. 

We have not yet determined under what conditions the plasma oscilla- 
tions can be considered to be in a frozen or non-frozen flow. Until this 
determination is made the time scale of coherence has some ambiguity 
associated with it, because our measured values include both the modifica- 
tion of the turbulence and the effect of the oscillations moving past the 
measuring point. It should be noted here that if the convection of the 
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oscillations past the probe is slow in comparison to the time it takes for 
an oscillation to be significantly modified, the time scales of coherence 1 
quoted will be very close to the true value. 

The measurements necessary to decide whether or not the plasma flow is 
frozen or continually being created and destroyed and those measurements 
needed to determine the convection velocity of the oscillation are being 
carried out at the present time. 

VI. Length Scales of Coherence and Strength of Turbulence 

The concept of a time scale of coherence was defined as a quantity 
which leads to an order of magnitude estimate of how long an oscillation 
remains in a specific form. Another quantity of interest is the length 
scale of coherence which determines over what spatial extent of the fluid 
the eddy or oscillation retains the same general form before it is signi- 
ficantly modified. 

The length scale of coherence, 1^ is defined in much the same way as 

the time scale except the definition of 1 utilizes the cross correlation 

function. The length scale of coherence is defined as the length over 

which T j x) at ^ = 0 drops to l/e of its maximum value^. Preliminary 

measurements taken in the radial direction have shown that 1 varies as much 

c 

as an order of magnitude in the frequency range 15 KHz - 100 KHz. 

*Y c an< 3- are quantities which involve the spatial and temporal extent 
of the turbulent oscillations. In order to completely describe the turbulent 
behavior of the plasma one needs a quantity which will indicate the strength 
or level of turbulence. The ratio of the a.c. fluctuation in charged particle 
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density n' to the measured d.c. charged particle density n Q is one such 
quantity. In our hollow cathode discharge n '/ n 0 varies from 0.1 % up to 
ICffo, depending upon the externally set operating conditions of the arc. 

Another term used in discussion of hydrodynamic turbulence is the 
"intensity of turbulence", which is defined as the ratio of the rms velo- 
city of the flow to the mean velocity 10 . 



where and 1 are the time and length scales of coherence respectively 

c c 

and v is the kinetic viscosity of the plasma. This quantity describes 
the plasma turbulent phenomena in terms of fluid turbulence terminology, 
time scale of coherence, power spectrum, length scale of coherence and an 
intensity of turbulence. Discussion of our results therefore requires 
consideration of the above parameters, and also whether we are dealing with 
collisional or collisionless plasma for comparison with theoretical predic- 
tions. Clarification of these points requires further work. 

Reports and Publications 

Dr. Noon reported on the progress of the work at NASA Lewis Center, 
Cleveland, in April 1970. A paper has been prepared for the November meeting 
of the Plasma Physics Division of the American Physical Society: "Correla- 

tion Measurements of Plasma Turbulence Spectra", G. Huchital, J. H. Noon, 
and W. C. Jennings. 
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